Dissipative maps at the chaos threshold: 
Numerical results for the single-site map 
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Abstract 

We numerically study, at the edge of chaos, the behaviour of the single-site map 
xt+\ = xt — xt/(xf + 7 2 ), where 7 is the map parameter. 
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In recent years, a number of studies has addressed the nonlinear dynamical 
systems at the chaos threshold [1-8]. At this special point (and also at some 
others such as period doublings, tangent bifurcations points where the Lya- 
punov exponent vanishes), it is shown that the proper entropy appears to be 
the nonextensive form [9] S q = k (l — Y^Li pf) /(<? — 1) with a proper q value 
(say q*) for each map under consideration. Up to now, three different methods 
are available for the determination of the q* values. First method comes from 
the analysis of the sensitivity to initial conditions where a general sensitivity 
function 

^t) = {l + (l-q)\ q t] 1/{1 - q) (qeR), (1) 



can be introduced. Here, X q is the generalized Lyapunov exponent and this 
form recovers the standard exponential one for q — 1, but it generically ex- 
hibits a power-law behaviour. In fact, eq.(l) corresponds to the power-law 
growth of the upper bound of the sensitivity function and this upper bound 
allows us to determine the proper q* value of the studied map. Second method 
is based on the geometrical aspects of the attractor at the edge of chaos and 
uses the scaling behaviour of the multifractal singularity spectrum f(a) as 
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where a m i n and a max being the most concentrated and most rarefied regions of 
the measure on the attractor respectively. It is obvious that calculating these 
values from the f(a) curve, one can easily estimate the proper q* value [3]. Fi- 
nally, third method, which was introduced very recently, deals with the entropy 
increase rates [5,6]. It basically based upon the (ensemble version) of gener- 
alized Kolmogorov-Sinai entropy K q = lim^oo limy^oo limjv^oo [^(t)] jt and 
it is conjectured that (i) a special q* value exists such that K q is finite for 
q = q*, vanishes for q > q* and diverges for q < q*; (ii) this q* value coincides 
with that coming from the other two methods explained above. 

Now we are prepared to introduce the main purpose of the present work. Here, 
we shall try to verify once again these conjectures using the abovementioned 
three methods for the single-site map 

x t+i = x t 9— — ? , (3) 
xf + 7 2 



where 7 is the map parameter [10]. The shape of the map is illustrated in 
Fig. 1. This map has a route to chaos by period doubling with Feigenbaum 
scaling but it exhibits an interesting property that the bifurcation cascade 
splits into two independent ones, which both obey Feigenbaum scenario. This 
property is the generic peculiarity for the single-site map in contrast to the 
logistic map [10]. On the other hand, these two independent cascades obey 
Feigenbaum scaling, which means that the Feigenbaum numbers of this map 
are the same as those of the logistic map, therefore, in that sense the single-site 
map shares the same universality class with the logistic map. Consequently, 
what we expect is that, at the chaos threshold, the proper q* value of this map 
would coincide with the value of the logistic case (namely, q* ~ 0.24). Our 
numerical results check and verify this expectation. Firstly, let us perform 
the first method to estimate the proper q* value of the single-site map at 
the onset of chaos (namely at 7 C = 0.43400483...). For the extremal point 
of the map (xq = 0.3213...), the sensitivity function exhibits a power-law 
divergence, £(t) oc t 1 ^ 1 ^^ as shown in Fig. 2. Measuring the upper bound 
slope 1/(1 — q), we determine the proper q value as q* = 0.24 . Now, let us 
proceed with the second method. In Fig. 3, the f(a) curve is given, where 
the monotonic behaviour with respect to the number of iterations is evident. 
Since the correct values of a m i n and a max could be obtained in principle when 
the number of iterations goes to infinity, the extrapolation shown in the inset 
of Fig. 3 is used to estimate the correct values of a min and a max . Using these 
values in the scaling relation Eq.(2), it is found that the proper q value is again 
q* = 0.24. (We also verify using the extrapolation technique that the fractal 
dimension df of this map coincides with that of the logistic map as expected). 
Finally, we shall use the third method for this map. To do this, we partition 
the phase space, x 6 (— 1/27, 1/27), into W equal cells, then we choose one 
of them and select N initial conditions all inside the chosen cell. As time 
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evolves, these N points spread within the phase space and finally this yields 
a set of probabilities. In the beginning of time, of course, S g (0) = 0, then 
it exhibits three successive regions. In the intermediate region, the entropy 
starts increasing and this is the region where the linear increase of the proper 
entropy is expected to emerge (for details see [5,6]). In Fig. 4, we plotted 
the time evolution of the entropy for three values of q. (At this point, it 
should be mentioned that the procedure of averaging over the efficient cells as 
described in [5,6] has also been performed here). In the intermediate region, it 
is seen that the linear increase of the entropy with time emerges for a special 
value of q. To support this picture, we fit the curve with the polynomial 
S q = A + Bt + Ct 2 in the intermediate region [ti,^]- Since the nonlinearity 
coefficient R = C(ti + t-^jB is a measure of the importance of the nonlinear 
term, it should vanish for a strickly linear fit. Therefore, inset of Fig. 4, we 
estimate the proper q value as q* = 0.234, which is the value when R is zero. 

Summing up, we study the behaviour of the single-site map at the chaos 
threshold. We numerically verify that the proper entropic index for this map 
is q* ~= 0.24. This constitutes another example of the one-dimensional dissi- 
pative maps which behave as q ^ 1 at the chaos threshold. 
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Figure Captions 

Figure 1 - The shape of the single-site map with 7 = 0.3 . 

Figure 2 - Log-log plot of the sensitivity function versus time. 

Figure 3 - Approximate multifractal singularity spectrum of the critical 
attractor of the single-site map for increasingly large number of iterations, going 
from N = 64 to 4096. Inset: a m i n , df and a max from the singularity spectra 
obtained from distinct number of iterations. Extrapolated values of them are 
also given. 

Figure 4 - Time evolution of S q (t) for three different values of q. W = 10 5 and 
N = 10 6 is used. The interval characterizing the intermediate region is [9,21]. 
Inset: The nonlinearity coefficient R versus q (see text for details). 
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